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RANGE  RESOLUTION  OF  TARGETS 


INTRODUCTION 

When  attempting  to  resolve  targets,  one  can  either  resolve  them  in  angle  or  range  or  a  combina¬ 
tion  of  both  coordinates.  The  problem  of  angular  resolution  has  received  considerable  attention,  but 
the  range  coordinate  is  a  much  more  powerful  discriminate  for  separating  targets,  for  example,  two  tar¬ 
gets  separated  by  1000  ft  are  separated  by  two  l-/xs  pulsewidths  but  by  only  1/20  of  a  2°  beamwidth  at  a 
range  of  100  nmi. 

The  probability  of  resolving  targets  is  not  only  a  function  of  target  separation  but  also  a  function 
of  their  signal  strengths  and  phase  differences.  However,  it  is  usually  assumed  that  if  targets  are 
detected,  they  can  be  resolved  if  they  are  separated  by  a  pulsewidth  [11.  In  practice  this  resolution 
capability  is  rarely  achieved.  For  instance,  a  commonly  used  resolution  algorithm  which  makes  no 
assumptions  about  the  shape  of  the  returned  pulse  does  not  attain  a  resolution  probability  of  0.9  until 
the  targets  are  separated  by  a  distance  greater  than  2.S  pulsewidths  [2],  Skolnik  [3]  conjectures  by 
drawing  an  analogy  with  angular  resolution  that  targets  within  0.8  pulsewidth  of  one  another  should  be 
resolved.  The  purpose  of  this  report  is  to  generate  and  evaluate  algorithms  for  achieving  resolution 
when  targets  are  within  a  pulsewidth  of  one  another. 

We  will  approach  the  resolution  problem  as  a  binary-hypothesis  problem  where  the  two  hypothesis 
are: 

•  Ht:  One  Target  Present 

•  H2:  Two  Targets  Present 

This  hypothesis  testing  approach  to  the  resolution  problem  was  first  considered  by  Helstrom  [4]  who 
investigated  the  maximum-likelihood  detector  for  the  detection  of  two  signals  of  known  separation  in 
time  and  frequency,  but  with  individual  phases  and  amplitudes  and  time  of  arrival  of  the  pair  unknown. 
Later,  Allen  [5]  considered  the  problem  of  resolving  targets  whose  positions  are  known  but  whose  pulse 
shape  may  be  distorted.  Our  approach  to  the  problem  will  be  the  generalized  likelihood  ratio  test  [6] 
where  the  unknown  parameters  are  replaced  by  their  maximum  likelihood  estimates.  Root  [7]  first  sug¬ 
gested  but  later  rejected  this  technique.  Although  the  approach  cannot  be  used  for  an  m- ary  hypothesis 
test  (H,:  exactly  /'  targets  present)  because  Hm  will  always  be  the  most  likely  hypothesis,  it  can  be  used 
for  the  binary-hypothesis  problem  yielding  a  desired  false-alarm  rate.  Then  the  m-array  hypothesis 
problem  can  be  attacked  by  first  performing  a  binary  hypothesis  as  to  whether  there  are  one  or  two  tar¬ 
gets  present.  If  one  decides  that  two  targets  are  present,  one  performs  another  hypothesis  test  deciding 
whether  two  or  three  targets  are  present.  One  continues  the  testing  until  the  null  hypothesis  is 
accepted. 

GENERALIZED  LIKELIHOOD  RATIO  TEST -NONCOHERENT  APPROACH 


We  now  formulate  the  resolution  problem  as  a  binary-hypothesis  test.  The  received  samples  are 


X,  -  n,  +  A  i 


sin  (a(R\  -  i)) 

- — - r —  co 

a(/?i  -  /) 

sin  (a(Ri  -  0)  , 


cos  9  +  A  2 


sin  (a(R2  -  0) 


a(R2-i) 
sin  (a(R2  -  /)) 


cos  (9  +  fi)  and 


r' '  ">■ +  A‘  sin  * +  Ai  „-r,—  sln  <(l + *’■ 
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where  X )  and  K,  are  the  inphase  and  quadrature  samples  for  the  ith  range  cell.  nx  and  nv  are  the  in¬ 
dependent  inphase  and  quadrature  Gaussian  noise  samples  with  mean  zero  and  variance  cr2,  A )  and  A2 
are  the  target  amplitudes,  Rt  and  R2  are  the  target  ranges  normalized  by  the  sampling  interval,  u  is  the 
phase  of  the  first  target,  f}  is  the  phase  difference  between  the  two  targets,  a  =  2. 7832/ A R,  and  A/?  is 
the  3-dB  pulsewidth  normalized  by  the  sampling  interval. 

We  have  assumed  that  the  targets  of  interest  are  point  targets.  This  assumption  will  not  affect  our 
results  as  long  as  the  extent  of  the  target  is  less  than  a  1/10  of  the  pulsewidth.  Hence,  the  results  of 
this  analysis  do  not  apply  to  high-range  resolution  radars.  The  range  resolution  problem  is  equivalent 
to  deciding  which  hypothesis  is  true: 

Hx  :  A  i  >  0  and  A  2  -  0 
or 

H2:  A  i  >  0  and  A2  -  0. 

If  {y4 ! ,^4 2,/? , ,/? 2. /3}  are  are  known,  the  optimal  detector  (in  the  Neyman  Pearson  sense)  given  the  /th 
range  sample  is 

A/  “  1 

27  Jo  P^.YMWdO 

where  P(Xh  T/lH;)  is  the  joint  density  of  Xt  and  Y{  given  that  hypothesis  H;  is  true.  Appendix  A 
shows  that  the  likelihood  ratio  A ,  reduces  to 

-Z,j/2<72  . 
e  /c 

A, - 

e  /  o 

where 

2}  -  A l  W}  +  2  AxA2  W,U,  cos  (3  +  A l  U,2, 

sin  (a(R  j  -  /)) 

a(/?i  -  i ) 

wi  “ 

0 


sin  (a(R 2  —  /)) 
a(R2  -  /) 

U,  - 

0 


,,  -  (jr;  +  Y2)'*2, 

and  /0(  )  is  the  modified  Bessel  function  of  the  first  kind.  Hence  the  likelihood  ratio  is  just  the  ratio  of 
two  Ricean  densities;  under  H2  there  are  two  signals  present  and  the  resultant  signal  is  given  by  Z,; 
under  H!  there  is  only  one  signal  present  with  amplitude  Wt  A,. 


|a(/?2  -  /) |  <  7T 
|o(Rj  —  f )  1  >  7T, 


la(/?i  -  i ) |  <  7 t 
|a(Rt  -  t)|  >  tr. 
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Given  the  envelope  detected  samples  rt,  i  —  1,  ....  m,  the  optimal  test  statistic  is 

A  —  fj  A(, 

t-i 

or  equivalently 


A'  -  log  A,. 

Since  in  the  problem  of  interest,  are  unknown,  we  replace  them  with  their  max¬ 

imum  likelihood  estimates.  Specifically,  we  find  the  values  of  the  parameters  which  minimize 

£  -  Y  ( X,-  JV/B  -  U,C)2  +  £  O',  -  W,D  -  U,E)2.  (2) 

/-i 

where  Z?  -  ^  cos  0,  C  - /f2  cos  (0  + /3),  £>  - sin  0,  and  E  -  A2  sin  (0  +  /3).  For  any  given 
value  of  Ri  and  R2,  ^  and  Ui  are  known  and  the  maximum  likelihood  estimates  of  B,  C,  D ,  and  E 
can  be  found  simply  by  solving  the  equations 

dL  SL  _  BL  6Z.  _  0 
a  B~  3C  "  dD~  be” 

The  solutions  for  the  one  target  case  (that  is,  assuming  H,  is  true)  are 

C  -  £  -  0, 

B  -  I  W,XJZ  W}.  (3) 

D-i  m T,/E  w?, 

where  all  sums  run  from  /  -  1  to  m.  The  solutions  for  the  two-target  case  (that  is,  assuming  H2  is 
true)  are 

B  -  1(1  Z/,2)(E  fVjXj)  -  (E  Wt  i/,)(E  U,X,))/ A, 

C  -  [(I  W?)( I  (/,*,)  -  (E  t/2) (E  W,X,)V A, 

Z)  -  l(E  £/,2)  (E  WJt  -  (E  ft',  £/,)  (E  U,  T,)]/A, 

£  -  [(E  W^2) (E  Z/,y;)  -  (E  N'T/, ME  f^fyl/A, 

where 

A  -  (E  ft',2) (E  Z/,2)  -  (E  WtU,)2. 

A  direct  search  is  used  to  estimate  the  target  locations  Z?i  and  R2.  Appendix  B  outlines  this  procedure. 
The  threshold  T  is  set  so  that 

Pr{\'  >  r|H,}  -  PFA~  a, 

where  PFA  is  the  probability  of  false  alarm.  Since  testing  for  multiple  targets  is  only  performed  after  a 
target  is  detected,  a  fairly  high  false-alarm  rate  is  permissible.  In  this  report  we  have  selected  a  - 
0.01.  Because  of  this  high  false-alarm  rate,  Monte  Carlo  simulations  can  be  used  to  estimate  £.  One 
thousand  repetitions  were  run  for  the  case  where  the  sampling  rate  was  1.5  samples  per  pulsewidth,  m 
-  8  samples,  cr  «*  1,  Ax  —  14.14  (20  dB)  and  A2  —  0.  For  each  repetition  the  pulse  location  was 
located  randomly  (uniformly  distributed)  with  respect  to  the  samples.  The  estimated  value  for  the 
threshold  was  3.7.  Essentially,  the  same  threshold  was  obtained  when  A\  was  raised  to  30  dB.  To 
obtain  information  on  the  relationship  between  resolution  and  sampling  rate,  a  higher  sampling  rate  of 
3.0  samples  per  pulsewidth  will  also  be  considered.  For  the  parameters  m  ”  10  samples,  <t  =  1,  A\  ” 
14.14  (20  dB),  and  A2  “  0,  the  estimated  threshold  was  3.5 
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PROBABILITY  OF  RESOLUTION 

The  probability  of  resolving  two  20-dB  nonfluctuating  targets  at  a  false-alarm  rate  of  a  =  0.01  was 
found  by  simulation  and  Figs.  1  and  2  show  the  results  as  a  function  of  target  separation  for  /3  =  0°, 
45°,  90°,  135°,  and  180°.  Two  targets  are  resolvable  if  they  are  separated  by  1/4  to  3/4  pulsewidths' 
depending  on  their  phase  relationship.  In  general,  it  is  easier  to  separate  targets  if  they  are  out  of  phase 
(that  is,  (8  close  to  180°).  The  reversal  of  the  135°  and  180°  curves  is  caused  by  the  fact  that  the  reso¬ 
lution  test  between  one  and  two  targets  is  not  applied  until  after  a  target  is  detected.  When  the  targets 
are  180°  out  of  phase,  there  is  a  large  target  cancellation.  In  fact  the  probability  of  resolution  curve  is 
practically  identical  with  the  probability  of  detection  curve;  that  is,  when  two  20-dB  targets  are  180°  out 
of  phase,  they  can  be  resolved  as  soon  as  they  can  be  detected.  The  dashed  curve  in  Fig.  2  represents 
the  resolution  capability  of  deciding  between  one  or  two  targets  when  the  detection  criterion  is  not 
imposed.  As  expected  this  curve  shows  that  it  is  easier  to  resolve  targets  when  —  180°  than  when  fi 
-  135°. 


Fig.  1  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A R  —  I.S  samples  per  pulsewidth;  target  strengths— nonfluctuating, 
Ax  —  A  2  —  20  dB;  phase  differences  —  0",  45°,  90°,  135°,  and  180°. 


Fig.  2  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A  R  -  3.0  samples  per  pulsewidth;  target  strengths— nonfluctuating, 
Ax  —  A2  “  20  dB;  phase  differences  “  0°,  45°,  90°,  135°,  and  180°. 


NRL  REPORT  8728 


Next,  we  investigated  the  effects  of  unequal  target  strengths.  Figures  3  and  4  show  the  results  for 
Ax  -  20  dB  and  A2  “  30  dB.  Figures  5  and  6  show  the  results  for  At  -  20  dB  and  A2  -  13  dB. 
Comparing  Figs.  1,  3,  and  5,  and  Figs.  2,  4,  and  6,  there  appears  to  be  less  variation  with  phase 
difference  when  the  targets  have  unequal  target  strengths.  The  A\  —  20  dB  and  A2  —  30  dB  cases  fall 
in  the  middle  of  the  Ax  -  A2  “  20  dB  cases,  whereas  the  A,  —  20  dB  and  A2  =  13  dB  cases  fall  to  the 
right  (worse  performance)  of  the  At  -  A2  —  20  dB  cases. 

I  next  investigated  the  effects  of  sampling  rate.  In  addition  to  the  1.5  and  3.0  samples-per- 
pulsewidth  case,  I  investigated  1.0,  2.0  and  4.0  samples-per-pulsewidth.  Repeating  the  procedure  for 
estimating  the  threshold,  thresholds  of  4.3  (m  “  7),  3.1  (m  =  8)  and  3.5  (m  “  14)  were  found  for  a 
false-alarm  rate  of  a  “  0.01  for  the  1.0,  2.0,  and  4.0  samples-per-pulsewidth  cases.  The  probability  of 
resolution  curves  are  shown  in  Figs.  7,  8  and  9.  Table  1  gives  the  range  separation  necessary  to  achieve 
a  0.5-resolution  probability.  As  the  sampling  rate  is  increased  from  1  to  4  samples  per  pulsewidth  (a 
400%  increase),  the  corresponding  decrease  in  range  separation  averages  only  41%.  Therefore,  while 
increasing  the  sampling  rate  improves  resolution,  it  does  not  appear  that  it  is  worth  the  cost  of  faster 
A/D  converters  and  increased  memory  and  processing.  The  best  compromise  appears  to  be  either  1.5 
or  2.0  samples  per  pulsewidth. 


Fig.  3  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A  R  “  1.5  samples  per  pulsewidth;  target  strengths— nonfluctuating,  A ,  — 
20  dB  and  A2  ~  30  dB;  phase  differences  “  0°,  45°,  90°,  135°,  and  180°. 


Fig.  4  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A  R  —  3.0  samples  per  pulsewidth.  target  strengths— nonfluctuating.  At  - 
20  dB  and  A2  —  30  dB;  phase  differences  —  0°,  45°,  90°.  135°,  and  180° 
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-ig.  5  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro- 
redure;  sampling  rate  A R  —  1.5  samples  per  pulsewidlh,  target  strengths— nonfluctuating,  At  — 
!0  dB  and  A 2  —  13  dB,  phase  differences  —  0°,  45°,  90°,  135°,  and  180°. 


-ig.  6  —  Probability  of  resolution  as  a  function  ol  range  separation:  Noncoherent  likelihood  pro- 
redure,  sampling  rate  A R  —  3.0  samples  per  pulsewidth;  target  strengths— nonflucluating,  A,  - 
!0  dB  and  A2  “  13  dB;  phase  differences  -  0°,  45°,  90°,  135°,  and  180°. 


Fig  7  —  Probability  of  resolution  as  a  function  of  range  separation  Noncoherent  likelihood  pro¬ 
cedure.  sampling  rate  AR  -  10  samples  per  pulsewidth.  target  strengths—  nonfiuctuating. 
At  —  A j  —  20  dB,  phase  differences  —  0°,  45°,  90°.  135°.  and  180° 
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Fig  8  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  \R  —  2.0  samples  per  pulsewidth.  target  strengths— nonfluctuating. 
A  |  —  /f  j  “20  dB,  phase  differences  “  0°,  45°,  90°,  135°.  and  180° 


Fig  9  —  Probability  of  resolution  as  a  function  of  range  separation  Noncoherent  likelihood  pro¬ 
cedure.  sampling  rate  \R  =  4.0  samples  pe:  pulsewidth.  target  strengths—  nonfluctuating. 
A  |  -  .4;  ~  20  JB;  phase  differences  =  O'.  45".  90"  135°,  and  180' 


Table  1  —  Range  Separation  Required  to  achieve 
a  Probability  of  Resolution  =  0.5 


Phase 

Difference 

(deg) 

Number  of  Samples  per  Pulsewidth 

(1.0) 

(1.5) 

(2.0) 

B19 

KSH 

■ 

Hfl 

■ 

1 

1 

EH 

■  ■ 

0.56 

0.40 

0.30 

45 

0.66 

0.51 

0.42 

0.37 

0.36 

0 

0.68 

0.54 

0  4$_ 

0.42 

0.39 
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MULTIPLE  PULSES 

I  now  consider  the  case  of  multiple  pulses  on  target.  There  are  several  situations  that  one  can 
consider.  These  situations  depend  upon  whether  the  parameters  lAi,A2,Rl,R2,0,fi)  remain  constant  or 
change  between  pulses.  Attention  will  be  restricted  to  the  situation  when  pulse-to-pulse  frequency  agil¬ 
ity  is  used.  In  this  case  it  is  reasonable  to  assume  that  Rt  and  R2  remain  constant  for  all  pulses  but 
that  A i,  A 2,  0,  and  (3  are  independent  per  pulse.  If  Xtj  and  Yu  are  the  yth  time  samples  of  the  /th 
range  sample,  it  is  straightforward  to  show  that  Bj  and  C,  can  be  found  from  Eqs.  (3)  and  (4)  by 
replacing  X,  with  X{j  and  similarly  Dt  and  Ej  can  be  found  from  Eqs.  (3)  and  (4)  by  replacing  Y,  with 

The  total  likelihood  ratio  is  just  the  product  of  the  individual  likelihood  ratios.  The  thresholds  for 
a  =  0.01  and  a  20-dB  nonfluctuating  target  for  n  *=  1,2,  4,  and  8  pulses  were  found  by  simulation  and 
are  given  in  Table  2. 


Table  2  —  Threshold  Yielding  a  Probability 
of  False  Resolution  =  0.01 


Number 

of 

Pulses 

Sampling  Rate 

1.5 

( m  =  8) 

3.0 

(m  =  10) 

n  =  1 

3.7 

3.5 

n  =  2 

4.7 

5.0 

n  -  4 

6.0 

6.3 

n  —  8 

8.5 

10.0 

Figures  10  and  11  show  the  results  for  two  20-dB  Rayleigh  fluctuating  targets.  The  phase 
difference  was  uniformly  distributed  between  0  and  180°.  The  reason  for  the  poor  performance  of  the 
one-pulse  case  is  because  the  probability  of  detection  is  only  0.9  for  a  20-dB  Rayleigh  fluctuating  target. 
Tables  3  and  4  give  the  ranges  at  which  two  targets  can  be  resolved  at  0.5  and  0.9  probabilities  when 
multiple  pulses  are  available. 


Fig  10  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure,  sampling  rate  AR  -  i.5  samples  per  pulsewidth.  target  strengthstcmfluctuating.  .4 ,  -  .4, 
-  20  dB;  phase  difference  are  independent  and  uniformly  distributed  per  pulse,  number  ol  pulses 
K  -  1.  2.  4.  and  8 
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Fig  1 1  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure.  sampling  rate  i/f  -  3.0  samples  per  pulsewidth;  target  strengthslemfluctualing.  .4,  -  A: 
-  20  dB:  phase  difference  are  independent  and  uniformly  distributed  per  pulse,  number  o'  raises 
,V  —  1,  2,  4.  and  8. 


Table  3  —  Range  Separation  Required  to  Achieve 
a  Probability  of  Resolution  =  0.5 


Number  of 
Pulses 

Number  of  Samples  per  Pulsewidth 

(1.5) 

(3.0) 

1 

o7?7  1 

0.42 

2 

0.28 

0.22 

4 

0.15 

0.14 

8 

0.10 

0.09 

Table  4  —  Range  Separation  Required  to  Achieve 
a  Probability  of'  Resolution  =  0.9 


Number  of 
Pulses 

Number  of  Samples  per  Pulsewidth 

(1.5) 

(3.0) 

2 

4 

I 

8 

0.18 

0.13 

AD  HOC  NONCOHERENT  PROCEDURES 

The  generalized  likelihood  test  involves  a  6-parameter  minimization  implemented  by  a  two  dimen¬ 
sional  search  (see  Appendix  B)  and  consequently,  is  very  time  consuming.  Therefore,  I  investigated  a 
combination  of  2  adhoc  approaches.  The  first  approach  involves  the  generation  of  double  peaks  and  the 
second  involves  two  small  of  a  maximum  amplitude  for  a  fixed  number  of  samples  above  the  detection 
threshold. 

Figure  12  shows  an  example  of  a  double  peak.  To  have  a  double  peak  it  is  necessary  for  the 
amplitudes  to  increase  (first  peak),  decrease  (null),  and  increase  (second  peak)  again.  Furthermore, 
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Fig.  12  —  Example  of  a  double  peak 


both  peaks  must  be  above  the  detection  threshold.  There  is  no  condition  on  the  null  between  the  two 
peaks.  When  the  sampling  rate  is  1.5  samples  per  pulsewidth,  a  double  peak  is  a  good  indication  of  two 
targets.  For  instance,  in  10,000  repetitions  of  a  20-dB  target,  no  double  peaks  were  observed.  On  the 
other  hand,  if  the  sampling  rate  is  3.0  samples  per  pulsewidth,  a  double  peak  is  no  longer  a  good  indica¬ 
tion  of  two  targets.  For  instance,  in  10,000  repetitions  of  a  20-dB  target,  442  double  peaks  were 
observed.  Since  this  corresponds  to  a  false-alarm  rate  of  0.044,  the  procedure  was  modified.  A 
modified  double  peak  is  the  same  as  a  double  peak  except  that  there  is  an  additional  condition  on  the 
null  depth.  Specifically,  the  null  depth  A N  defined  by 

AN  =  Minimum  {First  Peak  —  Null  Value,  Second  Peak  —  Null  Value) 

must  be  greater  than  a  specified  value.  The  required  null  depth  thresholds  were  found  by  simulation 
and  Table  5  gives  the  results.  When  there  are  multiple  pulses,  two  targets  are  declared  if  any  of  the 
pulses  contains  a  modified  double  peak. 

A  double  peak  is  most  likely  to  occur  when  the  two  targets  are  180°  out  of  phase.  On  the  other 
hand,  if  the  targets  are  in  phase,  one  would  expect  a  small  peak  value  for  a  given  number  of  samples 
greater  than  the  detection  threshold.  Simulations  were  run  to  generate  histograms  of  peak  return  for  a 
given  number  of  samples  above  the  detection  threshold.  For  instance,  for  10,000  repetitions  of  a  20-dB 
Rayleigh  fluctuating  target,  1302  times  there  were  no  samples  above  the  detection  threshold.  1210 
times  there  was  one  sample  above  the  threshold,  4703  times  there  were  two  samples  above  the  thres¬ 
hold,  2783  times  there  were  three  samples  above  the  threshold,  and  one  time  four  samples  were  above 
the  threshold.  Figure  13  shows  the  amplitude  histogram  when  there  are  two  samples  above  the  thres¬ 
hold.  If  the  resolution  threshold  is  set  at  5.95  a  false-alarm  rate  approximately  equal  to  0.01  is 
obtained.  To  obtain  the  resolution  threshold  when  there  are  four  samples  above  the  detection  thres¬ 
hold,  the  simulation  was  repeated  with  a  30-dB  fluctuating  target.  When  there  are  multiple  pulses  on 
target,  the  pulses  are  integrated  noncoherently;  then  the  same  procedure  is  applied  to  the  noncoherent 
sums.  The  results  for  1.5  samples  and  3.0  samples  per  pulsewidth  are  given  in  Tables  6  and  7  respec¬ 
tively. 

Using  the  resolution  thresholds  of  Tables  6  and  7,  the  probability  of  resolution  curves  were  gen¬ 
erated  and  are  shown  in  Figs.  14  and  15.  Comparing  the  adhoc  results  with  the  likelihood  approach 
results  appearing  in  Figs.  10  and  11,  one  observes  that  the  likelihood  approach  is  superior;  in  fact,  the 
likelihood  approach  for  2  pulses  at  a  sampling  rate  of  1.5  samples  per  pulsewidth  is  approximately  the 
same  as  the  adhoc  approach  for  8  pulses  at  a  sampling  rate  of  3.0  samples  per  pulsewidth. 
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Table  5  —  Thresholds  for  Null  Depth 
Necessary  to  Achieve  a  False  Alarm  Probability 
of  a  —  0.01 


Null  Depth  Threshold  Values 

Number 

of 

Number  of  Samples  per  Pulsewidth 

Pulses 

(1.5) 

(3.0) 

1 

0 

1.25 

2 

0 

1.40 

4 

0 

1.60 

8 

0 

2.00 

AMPLITUDE  OF  LARGEST  RETURN 

Fig.  13  —  Amplitude  histogram  of  case  where  there  are  two  pulses  about 
detection  threshold:  .4,  “20  dB  and  A2  =  -  °°  dB 


Table  6  —  Resolution  Threshold  Yielding 
a  False-Alarm  Rate  of  0.01  when  the  Sampling  Rate 
_ is  1.5  Samples  per  Pulsewidth 


Number  of 
Pulses 

Number  of  Samples  Above  Threshold 

(2) 

(3) 

(4) 

(5  and  above) 

1 

5.95 

7.35 

33.65 

OO 

2 

9.82 

15.32 

51.42 

oo 

4 

19.48 

28.78 

87.88 

OO 

8 

50.55 

62.55 

100.05 

oo 

Table  7  —  Resolution  Threshold  Yielding  a  False-Alarm  Rate 
of  0.01  when  the  Sampling  Rate  is  3.0  Samples  per  Pulsewidth 


Number  of 
Pulses 

Number  of  Samples  Above  Threshold 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8  and  above) 

1 

■ 

WEEM 

8.85 

46.75 

oo 

2 

14.92 

24.32 

56.42 

OO 

4 

1 

SH 

26.88 

36.38 

OO 

8 

- -  - 

■HI 

55.25 

68.05 

200.05 

oo 

H 


NOI 


GERARD  V  TRUNK 


Fig.  14  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A R  —  1.5  samples  per  pulsewidth,  target  strengthsfemfluctualing,  At  —  A2 
**  20  dB,  phase  difference  are  independent  and  uniformly  distributed  per  pulse,  number  of  pulses 
N-  1,2,  4,  and  8. 


Fig.  15  —  Probability  of  resolution  as  a  function  of  range  separation:  Noncoherent  likelihood  pro¬ 
cedure;  sampling  rate  A/?  “  3.0  samples  per  pulsewidth;  target  strengths (emfluctuating,  At  -  A2 
—  20  dB;  phase  difference  are  independent  and  uniformly  distributed  per  pulse,  number  of  pulses 
N  -  1,  2,  4,  and  8. 


Table  8  —  Threshold  for 
Coherent  Likelihood  Ratio  with 
a  Sampling  Rate  of 
1.5  Samples  per  Pulsewidth 


Number  of  Pulses 

Threshold 

1 

12.4 

2 

18.5 

4 

26.7 

8 

39.0 
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GENERALIZED  LIKELIHOOD  RATIO  TEST-COHERENT  APPROACH 

It  all  the  signal  parameters  except  «  are  known,  the  optimal  test  is  given  by  the 

likelihood  ratio  (Eq.  (1))  involving  the  noncoherent  samples.  Since  the  signal  parameters  are  un¬ 
known,  we  estimated  them  and  substituted  the  estimates  into  the  noncoherent  likelihood  ratio.  How¬ 
ever,  since  the  likelihood  ratio  of  noncoherent  samples  now  contains  estimates  instead  of  the  true 
parameters,  it  is  no  longer  optimal,  and  it  is  possible  that  the  likelihood  ratio  based  on  coherent  samples 
will  yield  better  results.  Consequently,  we  now  consider  the  coherent  likelihood  ratio. 

The  coherent  likelihood  ratio  for  n  pulses  is 


m 

U 

27exp 

-y  (X0  -  W,BJ  -  +  (Y,j  -  W.D,  -  U.Ef) 

m 

Jl 

2tt"  6XP 

«*„-  WiBf  +  iY^-  W^j)1 1 

Equation  (4)  gives  the  maximum  likelihood  estimates  for  the  two  target  signal  parameters  ( Bj.Cj.Dj.Ej ) 
appearing  in  the  numerator,  and  Eq.  (3)  gives  the  maximum  likelihood  estimates  for  the  single  target 
signal  parameters  [Bj.Dj]  appearing  in  the  denominator.  Taking  the  log  of  the  coherent  likelihood 
yields  as  a  test  statistic,  the  square  error  residue  (Eq.  (2))  evaluated  for  one  target  minus  the  square 
error  residue  (Eq.  (2))  evaluated  for  two  targets 

log  A  =  Min  [L]  -  Min  {/.}.  (5) 

R^.Bj.Dj  Rl,R1.Bj.CJ.DJ.Ej 

We  now  proceed  as  we  did  for  the  noncoherent  likelihood.  We  first  calculate  the  thresholds  for 
multiple  pulses  by  using  simulation  techniques.  The  thresholds  for  n  =  1,2,  4,  or  8  pulses  are  given 
in  Table  8.  Next,  again  using  simulation,  Figs.  16-19  show  the  resolution  curves  for  A,  -  A2  =  20  dB. 
A i  =  20  dB  and  A2  =  30  dB,  A\  =  20  dB  and  A2  =  13  dB,  and  multiple  pulses.  Comparing  the 
coherent  results  with  the  noncoherent  results  of  Figs.  1,  3,  5,  and  10,  one  concludes  that  the  coherent 
results  are  uniformly  better  and  are  much  better  when  the  targets  have  different  signal  strengths. 


Fig.  16  —  Probability  of  resolution  as  a  function  of  range  separation:  Coherent  likelihood  pro¬ 
cedure;  sampling  rate  if?  —  1.5  samples  per  pulsewidth;  target  strengths— nonfluctuating, 
A]  -  A i  “  20  dB;  phase  differences  “  0°,  45°,  90°,  135°,  and  180°. 
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AD  HOC  COHERENT  PROCEDURE 

Although  the  previously  discussed  ad  hoc  noncoherent  procedures  are  not  very  applicable  to  the 
coherent  data,  the  coherent  likelihood  ratio  suggests  another  ad  hoc  approach:  instead  of  using  the 
difference  in  square  errors  given  by  Eq.  (5),  use  the  square  error  associated  with  fitting  one  target  to 
the  data.  If  only  one  target  is  present,  the  residue  should  only  be  noise  and  hence  should  be  small.  If 
two  or  more  targets  are  present,  the  residue  will  contain  signal  and  should  be  large.  This  statistic  only 
requires  a  one-dimensional  search  and  can  be  accomplished  very  quickly.  In  fact  since  the  initialization 
procedure  given  in  Appendix  B  is  fairly  accurate,  it  is  unclear  whether  or  not  the  one-dimensional 
search  is  required. 

Using  simulation  techniques,  the  appropriate  thresholds  were  calculated  and  are  given  in  Table  9. 
Next,  again  using  simulation,  the  resolution  curves  for  At  —  A  2  =  20  dB,  A\  =  20  dB  and  A2  “  30 
dB,  A\  =  20  dB  and  A2  “  13  dB,  and  multiple  pulses  were  calculated  and  are  shown  in  Figs.  20  to  23. 
Comparing  the  ad  hoc  coherent  results  with  the  coherent  likelihood  results  shown  in  Figs.  16  to  19,  one 
concludes  that  the  coherent  likelihood  results  are  only  slightly  better.  Consequently,  because  of  the 
reduced  calculations  required,  I  would  recommend  the  ad  hoc  approach 


Table  9  —  Thresholds  for 
Fitting  Single  Target  to  the  Data 
with  a  Sampling  Rate  of 
1.5  Samples  per  Pulsewidth 


Number  of  Pulses 

Threshold 

h  1 

28.5 

2 

44.7 

4 

73.7 

8 

148.4 

Fig  20  —  Probability  of  resolution  as  a  function  of  range  separation:  Coherent  likelihood  pro¬ 
cedure;  sampling  rate  AR  —  1.5  samples  per  pulsewidth.  target  strengths— nonfluctuating, 
/t,  -  A 2  -  20  dB;  phase  differences  -  0°.  45°.  90°.  135°.  and  180° 
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O.l  0.2  0.3  0.4  0.5  0.6  0.7 

TARGET  SEPARATION  IN  PULSEWIDTHS 

?ig.  21  —  Probability  of  resolution  as  a  function  of  range  separation:  Coherent  likelihood  pro- 
edure;  sampling  rate  A R  —  l.S  samples  per  pulsewidlh;  target  strengths— nonfluctuating,  A2  — 
10  dB  and  A2  —  30  dB;  phase  differences  —  0°,  45°,  90°,  135°,  and  180°. 


O.l  0  2  03  0.4  0.5  0.6  0.7  0 

TARGET  SEPARATION  IN  PULSEWIDTHS 

■ig.  22  —  Probability  of  resolution  as  a  function  of  range  separation:  Coherent  likelihood  pro- 
edure;  sampling  rate  A R  —  1.5  samples  per  pulsewidth;  target  strengths— nonfluctuating.  A,  — 
!0  dB  and  A2  —  13  dB;  phase  differences  —  0°,  45°,  90°,  135°,  and  180°. 


0  4  0.5  0  6  0,7  0 

TARGET  SEPARATION  IN  PULSEWIDTHS 

Fig,  23  —  Probability  of  resolution  as  a  function  of  range  separation:  Coherent  likelihood  pro¬ 
cedure,  sampling  rate  A#  —  1.5  samples  per  pulsewidth;  target  strengthsfemfluctuating,  4,  «  A2 
—  20  dB;  phase  difference  are  independent  and  uniformly  distributed  per  pulse,  number  of  pulses 
N-  1,2,  4,  and  8. 
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SUMMARY 

The  problem  of  resolving  targets  in  range  has  been  formulated  as  a  hypothesis-testing  problem 
where  the  unknown  parameters  (target’s  location,  amplitude,  and  phase)  are  replaced  by  their  least- 
square  estimates.  When  the  likelihood  ratio  test  using  either  coherent  or  noncoherent  samples  is 
applied  to  the  range  samples  surrounding  a  detection,  good  resolution  results  are  obtained.  The  proba¬ 
bility  of  resolution  is  a  function  of  the  sampling  rate,  target  amplitudes,  target  separation,  and  target 
phase  difference.  Using  a  sampling  rate  of  1.5  samples  per  pulsewidth,  two  20-dB  nonfluctuating  targets 
can  be  resolved  at  a  resolution  probability  of  0.9  at  separations  varying  between  1/4  and  3/4  of 
pulsewidth  depending  on  the  relative  phase  difference  between  the  two  targets.  This  result  can  be 
improved  further  by  processing  multiple  pulses.  The  coherent  likelihood  is  preferred  because  it  yields 
better  results  when  the  two  targets  have  unequal  strengths. 

The  likelihood  approach  was  compared  to  an  easily  implemented  ad  hoc  approach.  The  ad  hoc 
approach  declared  the  presence  of  two  targets  if  either  of  two  conditions  were  met:  (1)  the  return  signal 
was  double  peaked  with  both  peaks  above  the  detection  threshold  or  (2)  the  peak  amplitude  was  too 
small  in  relationship  to  the  number  of  samples  above  the  detection  threshold. 

Unfortunately,  while  this  ad  hoc  approach  is  easily  implemented,  the  likelihood  approach  has 
superior  performance.  A  slightly  more  complicated  approach  involves  fitting  a  pulse  shape  to  the  data 
and  comparing  the  residue  error  to  an  appropriate  threshold.  This  procedure  is  only  slightly  less  accu¬ 
rate  than  the  coherent  likelihood  ratio  and  hence  is  the  recommended  procedure. 
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The  likelihood  ratio  is 


<p(x,y  |H2)>, 


where  the  subscript  for  the  ith  sample  has  been  dropped  and  the  notation  <  •>«  indicates  that  the  den¬ 
sity  is  averaged  over  0.  We  now  calculate  <p(x,.y  |H2)>«  and  obtain  </>(x,y|Hi)>8  from  it  by  setting 
the  amplitude  A2  of  the  second  signal  to  zero. 

The  inphase  and  quadrature  signals  x  and  y  can  be  written  as 

x  “  nx  +  A  i  W  cos  9  +  A  2  U  cos  (0  +  / 3) 
y  -  riy  +  A !  W  sin  9  +  A2  U  sin  (9  +  ft) 


where 


fV 


sin(a(Xi-0)_\a(Ri_i)\<ir 


a(R  i  -  i) 


|a(/?j  -  i)|  >  w 


U  - 


sinUtt.-ft)  |a(.  _„U, 


a(R  2-  i) 


®  —  /)  l  >  7T 

Since  x and  y  are  independent,  </»(x,y  |H2)>9  can  be  written  as 

<p(x,y \H2)>9  —  — exp  {  -  [x  -  A\  W  cos  9  -  A2U  cos  (0  +  /3)l2/2or: 

2n<r 

-  Ly  “  A\  W  sin  0  -  A2U  sin  (0  +  /3)]2/2<rJ}>9. 

Performing  the  squaring  operation  and  combining  terms  yields 

1 


<p(x,y|H2)>#  - 


2ir<r2 


exp  {  -  I x2  +  y2  +  AfW1  +  2 A\A2WU  cos  0 


+  AjU2]/2o-2)  <  exp  (-4-  lA,  Wx  cos  0 
<r 

+  A\ Wy  sin  0  +  A2Ux  cos  (0  +  0)  +  A2Uy  sin  (0  +  /3)]}># 
Letting  x  —  r  cos  <t>  and  y  —  r  sin  ^  yields 

<p(r,$|H2)>9  -  C  exp{-^-  lAiW  cos  (0  -  <A)  +  A2U  cos  (0  +  / 3  -  0)]}>8 
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where 

C  -  exp{  -  [r2  +  AfW2  +  2AXA2WV  cos/3  +  A^U2]/2ct2\ 

2ncr 1 

Since  the  integral  over  9  is  from  0  to  2jt,  the  value  of  the  integral  is  independent  of  0  and  thus  we  can 
set  0-0: 

<p(r,<t> |H2)>«  -  C<exp  { rlAt  W  cos  9  +  A2U  cos  (9  +  >3) ]/cr2} > s 
Expanding  cos  (9  +  0)  and  rearranging  yields 

<p(r,0|H2)>o  -  C<exp  {r[A\  W  +  A2U  cos  /3)cos  9  -  A2U  sin  /3  sin  0]/<r2}># 

Letting  {A\  W  +  A2V  cos  /3)  -  Z  cos  a  and  ~A2U  sin  0  -  Z  sin  a  yields 

</>(r,0|H2)>8  -  C<exp{rZ  cos  (0  -  a)/a2)>9 

The  preceding  expression  is  independent  of  a.  Setting  a  —  0  yields  an  integral  which  is  a  modified 
Bessel  function: 

<p(r,0|H2)>s  -  - r -r  exp  { — lr2  4-  Z2l/2cr2}  /0  ( rz/tr 2) 

2rt<rl 

where  Z 2  -  AftV2  +  2 A,A2WU  cos/3  +  A22U2. 

Integrating  over  0  and  substituting  into  the  likelihood  ratio  yields 

exp(—  Z2/ 2a2)  I0  irZ/a2) 
exp{-  A 2  H/2/2<r2}  /0  (rA\Whr2) 


Appendix  B 

SEARCH  PROCEDURE 


The  square  error  is  minimized  by  using  a  direct  search*  for  the  target  locations.  The  initial  start¬ 
ing  location  for  the  first  target  is  found  by  interpolating  around  the  largest  return.  To  do  the  necessary 
interpolation,  we  assume  that  the  pulse  shape  can  be  approximated  by  a  Gaussian-shaped  pattern 

A(R)  -  e-a(R-Rr)\ 

where  A  (R )  is  the  pulse  gain  at  range  R  when  a  point  target  is  located  at  range  Rr.  The  appropriate 
normalizing  constant  a  can  be  found  from 

0.5  -  le~aUW2)2]2, 


yielding 


2  ln(.S)  1.386 

R  IdB  R  3  dB 


Now,  assume  that  the  target  is  detected  and  the  corresponding  ranges  to  the  largest  return  and  next 
largest  adjacent  return  are  /?!  and  Rx  +  A/?  where  A/?  is  the  sampling  rate.  Then,  the  corresponding 
amplitude  returns  are  given  by 

A\  —  e  ' 


and 


*2 


Solving  these  equations  for  the  target  range  yields 

*r-*i  +  ¥+  why 108  w-*')' 

However,  the  above  estimate  need  not  lie  between  R\  and  R\  +  A R.  Consequently,  if 


and  if 


Rj  ^  R\  set  Rf  m 


Rt>  R}  +  A/?,  set  RT  -  /?!  +  A R. 


To  estimate  the  starting  location  for  two  targets  we  first  use  the  previous  procedure  to  find  the 
location  of  the  first  target.  We  then  noncoherently  subtract  out  the  first  target  and  apply  the  previous 
procedure  to  the  residue. 


•R.  Hooke  and  T.  Jeeves,  "Direct  search  solution  of  numerical  and  statistical  problems."  J.  ACM,  8.  212  (1961). 


